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Abstract 

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed 
kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. 
These times are then correlated to the reactor conditions using an exponential fit Simple first order reaction expressions are 
then used to find the conversion in the reactor. The method uses a two time step kinetic scheme. The first time averaged step 
is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of 
initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations 
(>lxlO -20 moles/cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water 
mole concentrations, pressure and temperature ( T 4 ). The simple correlations are then compared to the turbulent mixing times 
to determine the limiting properties of the reaction. 

The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic 
scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is 
regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H 2 and 
NO x are obtained for H 2 /Air fuel and for H 2 /Q 2 . 

A similar correlation is also developed using data from NASA’s Chemical Equilibrium Applications (CEA) code to 
determine the equilibrium temperature (T 4 ) as a function of overall fuel/air ratio, pressure and initial temperature (F 3 ). High 
values of the regression coefficient R 2 are obtained. 


Introduction 

Hydrogen is the fuel of the future. Iceland is converting their economy to hydrogen and several countries (United States, 
Germany, Australia, Mexico, China, Canada, Japan, and Russia) are strongly interested in the technology because of global 
warming concerns from CO 2 and the unlimited supply of hydrogen from water conversion. It is important to have a simple, 
quick method of computing the chemical rates for prediction of times of conversion and for the design of combustors. 
Accurate chemical kinetic mechanisms exist for hydrogen reaction, but they involve as many as 21 chemical steps with the 
production of intermediate species like OH, H0 2 , H and O atoms with H 2 02 as an intermediary. We have used the H 2 
mechanism of S.L. Hwang, et al. (ref. 1) as the basis of this report. No changes in the hydrogen mechanism were made in 
switching from H 2 /Air to H 2 /0 2 . We have taken a fourty-two step mechanism for the production of nitrogen oxides from 
GRI-Mech 2.1 without the carbon reactions, used in reference 2, and correlated the formation rate of nitrogen oxides in the 
presence of hydrogen by a similar two time step method. Any detailed mechanism could have been used to compute the 
chemical kinetic time constants. 

The hydrogen reactions are about seven times as fast as hydrocarbon fuels and the reactions require small time steps to 
integrate accurately. We have used the NASA Glenn GLSENS kinetics code (ref. 3), to integrate the progress of the reaction 
accurately. The classical methods were used to compute over two thousand cases at constant T and P. Then the results were 
correlated with a simple two time step approach to predict the chemical kinetic time constants. The derived method involves 
only the principal species H 2 , and H 2 0 and NO* for nitrogen oxides, species which are predicted routinely by computer 
codes. We do not require tracking the intermediate free radicals, as the kinetic times in step one before the radicals become 
large are practically constant and in step two are related to the instantaneous H 2 , 0 2 , and H 2 0 concentrations . 

This report is a summary of a more detailed summary report (ref. 4). The two time step approach was used successfully 
by Molnar and Marek in reference 2 for Jet- A and methane. The two time step approach is used because the free radicals OH, 
etc., catalyze the reactions and act as important intermediates. The principle reaction for water production is OH through the 
reaction: 


H 2 + OH = H 2 0 + H 

We think that the more water that is present in the reactor then the more OH would be present accelerating all of the 
reactions. So H 2 0 can be used as an indicator of the OH present. 

In order to obtain the chemical source term, a comparison can be made of the turbulent mixing time t m or kJAz with k 
being the turbulent kinetic energy, e is the dissipation rate and the chemical kinetic time x Ct and the longest time is used as the 
overall time (ref. 5). This may also be represented by the following relationship: 
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With the approach derived here, a simple direct comparison can be made between the mixing and chemical kinetic times and 
the minimum rate used for all points and at each iteration of a CFD code. 

Simple Model Equations 

Only the following equations can be used to model the chemical system. 


X Fuel 

H2 + V2 O2 ^ H2O 

%o 

V 2 N 2 + l /2 0 2 ► NO 

The following first order reaction is used to represent the rate of fuel burning. We compute the net reaction time so we 
do not worry about the forward and reverse processes, only the net process. We assume that the fuel goes to zero at the end of 
the reaction and that the NO x increases linearly and it is far from equilibrium. (In this report, t and x are given in 
milliseconds, while concentrations are given in gmoles/cc): 
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The chemical time X H 2 can be computed from the concentration of H 2 and its net rate of production. 

The nitrogen oxide formation rate, a species important for combustor emissions, was modeled as a simple zero order 
expression. 
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where x NO has units of 

x gmol 

H 2 /Air Initial conditions for chemical kinetic 

Temperature = 1000, 1500, 2000, 2500 K 

Pressure = 1, 10, 20, 30, 40 atm 

Lean equivalence ratios = 0.3 to 1.0 increments of 0.1 

Rich equivalence ratios = 1.05, 1.1 to 2.0 increments of 0.1 

H 2 /O 2 Initial conditions for chemical kinetic times for advanced rocket combustors (ref. 6) 

Temperature = 1000, to 4000, increments of 500, then 5000 K 
Pressure = 0.5, 1.0, 10, 50, 100, 200, 300, 400 atm 
Lean equivalence ratio = 0.1 to 1.0 increments of 0.1 
Rich equivalence ratio = 1.05, 1.1, 1.2 to 3.0 increments of 0.2 

Time steps for Step 2 were the same as Step 1 with integration to 3 milliseconds. The output tabulated at least 20 points 
per case. For the approximately 2000 cases, the total number of points correlated could exceed 40000 per regression. 
Excel can handle up to 65000 points. 

The kinetic mechanism of Hwang, S.M., et al. (ref. 1), next page, was used for the full range of conditions without 
modification. 

Even though the mechanism was tested to a temperature of 3000 K and 4 atmospheres, we felt that the mechanism would 
extrapolate and give reasonable results for design because nothing else was available at the high temperature and high 
pressure conditions. The differences between the H 2 /Air and the H 2 /0 2 is the absence of nitrogen in the H 2 /0 2 set, and in 
H 2 /0 2 the maximum temperatures were twice that of H 2 /Air and the maximum pressures were an order of magnitude higher. 
We computed and correlated the chemical kinetic times over the full range of conditions. 
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H 2 Kinetic Scheme for Step l(Average) and Step 2 (Instantaneous) Methods 




CC S.M. Hwang, et al. (ref. 1) H 2 Combustion Over a Range of Conditions 
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Some reactions are bimolecular and some are triinolecular expressions so the kinetics code 

| GLSENS was used. 







Correlation Results 

Equilibrium Lean and Rich H^Air Correlations 

The H 2 equilibrium correlation for T 4 was correlated for both the lean and rich cases (ref. 7). The exponential function 
did not correlate the complete data over the total equivalence ratio from 0.2 to 2.0, because of the peak in the values at the 
stoichiometric conditions. We chose to correlate the values from lean (<|> = 0.2 to 1.) and then rich from >1.05 to 2.0). A 
similar approach was used for the chemical kinetic times. The T 4 equations were correlated over the complete temperature, 
pressure and equivalence ratio O: The stoichiometric d> was included on the lean side. 
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Although these equations seem complex, they are simple to program and use. Calculations can even be done on a hand 
calculator. 



A 

a 

b 

c 

R-squared 

Lean 7 4 
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0.00419 

0.261 

0.432 

0.936 

Rich T 4 

1217.0 

0.00727 

0.120 

-0.243 

0.962 


R-square values have been included to demonstrate the strength of the correlation. R-squared is a measure of the error in 
the model; an R-squared value of one is ideal. The good R 2 values obtained are typical of all the variables. 


Chemical Kinetic Times for H^Air 

Step One Equations for H^Air Chemical Kinetic Time Correlations 


The following form of equation was used for the H 2 and NO x step one correlations: 

eT 

T _ 

where x is the chemical kinetic time in milliseconds, P is pressure in atm, #a is the initial or overall mass fuel air ratio as in 
reference 2, and T is the temperature in Kelvin. The correlation is switched to Step Two when the water molar concentration 
is greater than lxlCT 20 moles/cc. This value was chosen arbitrarily as a simple switch. Step two can not be used with small 
concentrations of water because if the value of the water concentration were zero, the entire correlation time would go to zero 
or predict an infinite rate. 


x = A(p) a (of exp 


Step One H 2 /Air 
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Step Two Equations for H/Air 


The following form of equation was used for both the x m step two fuel and NO* correlations except for different 
constants. This form of the equation produced the best fit: 


T = A{pf (cH 2 ) b (cH 2 0) c exp 


£ 

T 


where P is pressure in atm, cH 2 is the instantaneous molar concentration of H 2 and cH 2 0 is the instantaneous water 
concentration, moles/cc, and T is the mixture temperature in Kelvin. We have correlated to only the major species hoping that 
cH 2 0 will track the minor species (OH, H, O, etc.) to allow good overall correlation and easy to use equations. 


Step Two Ha/Air Chemical Kinetic Time Correlations 
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H2/O2 Equilibrium Correlations 

The H 2 /0 2 presentation closely follows the H 2 /Air presentation. 
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Equilibrium Lean and Rich H2/O2 Correlations 
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Again a good correlation is obtained. 


Chemical Kinetic Times for H2/O2 


Step One Equations for HJ0 2 


The lean and rich step one H 2 /0 2 chemical kinetic time correlations are of the following form: 


T = 


A(P) a (Phi) b ex p 


D 

T 


where P is pressure in atm, phi is the initial or overall equivalence ratio and T is the temperature in Kelvin. 


Step One Hj/O? Chemical Kinetic Time Correlations 
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Step Two Equations for H2/O2 


The following form of equation was used for all lean and rich H 2 /0 2 step two correlations: 
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Step Two H2/O2 Chemical Kinetic Time Correlations 
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Comparison of the Model With the Full Mechanism Kinetics 

The predictions of the full mechanism (H 2 , H, O, OH and H 2 0) are plotted in figure 1, and then our model predictions are 
superimposed on the calculations. The model predictions are close to the full mechanism predictions. The model used a very 
simple mass balance for the water concentration. We have not included the intermediate species in our predictions, but we 
correlate the chemical kinetic times to the overall species concentrations The H 2 0 model concentration is higher than the full 
mechanism, because no intermediate H species exist in the model. The ignition delay time (3.E-^4 milliseconds) is closely 
predicted. As expected, OH, H and O increase at almost the same rate as H 2 0, so H 2 0 is a good indicator of the reaction rate 
Since we use the model predictions of concentration to predict the chemical kinetic times, we have correlated the full 
mechanism chemical kinetic times against the model concentrations of H 2 and H 2 0. The chemical kinetic time for this 
condition is about one microsecond so the integration requires only about 3000 steps. 

Reduced Model used 

These simple equations in an Excel spreadsheet. 
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cH 2n = cH 2 „.i*INEFF 
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cH 2 0 = (cH 2 °-cH 2 ) 


cN° x/l = cNO x „_, 


tn—tn - 1 


Used backward differencing At = 1 .e-6 for this comparison. 

Figure 1. H2/Air Concentration (Full Mechanism) 
Compared to Reduced Model vs Time ms 
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Conclusions 

A simplified kinetic scheme for H 2 /Air and H 2 /0 2 fuels resulted in a two time step correlation that calculates chemical 
kinetic times for fuel and NO x . These models are simple enough to use in spreadsheet calculations and they are very accurate. 
These chemical kinetic time equations can then be used in a numerical combustor code to compare the chemical kinetic time 
with the turbulent mixing time at every grid point and for every iteration step. Strong step one and step two H 2 correlations 
were developed (R 2 > 0.9) for all relations, both the H 2 /Air and H 2 /0 2 correlations. However, because we are trying to 
correlate so many values over a wide range of conditions, error exists. These model equations are believed to be extremely 
useful in the comparison of kinetic reaction and turbulent mixing times and in the computation of kinetic rate results. In our 
model we have only two equations and four major species, therefore the step size can be large and the rates are no longer 
coupled. 
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